Consequences of Heavy Quark Symmetries for Hadronic Molecules 
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Among the newly observed structures in the heavy quarkonium mass region, some are proposed 
to be hadronic molecules. We investigate the consequences of heavy quark flavor symmetry on these 
heavy meson hadronic molecules. The symmetry allows us to predict new hadronic molecules on one 
hand, and test the hadronic molecular assumption of the observed structures on the other hand. We 
explore the consequences of the flavor symmetry assuming the X(3872) and ^(,(10610) as a isoscalar 
DD* and isovector BB* hadronic molecule, respectively. A series of hadronic molecules composed 
of heavy mesons are predicted. In particular, there is an isoscalar 1+"*" BB* bound state with a mass 
of 10580 MeV which may be searched for in the T(15', 2S')7r^7r~7r'' mass distribution; the isovector 
charmonium partners of the ^(,(10610) and the Zi,(10650) are also predicted, one of which probably 
corresponds to the very recently observed Zc(3900) resonance by the BESIII Collaboration. 

PACS numbers: 12.39.Hg, 14.40.Rt, 14.40.Pq, 13.75.Lb, 03.65.Ge 



Due to ongoing experimental efforts, a series of new 
heavy quarkonium states, called XYZ states, have been 
observed in the last decade. Many of them arc expected 
to be of exotic nature, for a comprehensive review wc re- 
fer to Ref. [H . Due to the proximity of the masses to cer- 
tain hadronic thresholds, some of the XYZ states have 
been proposed to be hadronic molecules, i.e. states that 
are generated by the interaction between two or more 
hadrons (they are bound states if they are below the 
threshold and in the first Riemann sheet, or virtual states 
and resonances if they are in the second Riemann sheet 
of the scattering amplitude). For instance, the famous 
X(3872) discovered by the Belle Collaboration 0] and 
confirmed by many other experiments was proposed to 
be a DD* (the charge conjugated particles are implicitly 
included here and in the following) bound state with an 
extremely small binding energy Q or a virtual state Q ; 
the isovector ^^'"(10610) and Z^'^10650) resonances 
reported by the Belle Collaboration 0, Q have been con- 
sidered to be BB* and B*B* hadronic molecules, re- 
spectively 0; the y(4660) observed in the V^'tttt mass 
distribution by the Belle Collaboration and confirmed 
by BaBar Q is possibly a ■(/)'/o(980) bound state [l^l; the 
X(4260) state observed by the BaBar Collaboration [llj 
has been suggested to be a DDi molecule [H, [13] . Other 
models also exist for these states. Thus, in order to 
understand these states and furthermore their binding 
mechanisms, it is important to find out methods which 
can distinguish the different scenarios. Decay patterns 
are often used for this purpose, here we will pursue a 
different approach using heavy fiavor symmetry. 



Without developing complicated dynamical models, 
symmetries are often helpful in describing certain as- 
pects of various physical systems. For a system in- 
volving a heavy quark whose mass mq is much larger 
than Aqcd, flavor and spin symmetries arise by sending 
toq to infinity (for reviews of heavy quark symmetries, 
see Refs. [3, 111). Due to spin symmetry, both heavy 
mesons and heavy quarkonia form spin multiplets, e.g. 
the {D,D*} and {"qcJ/il}}. The masses are degenerate 
in the heavy quark limit, and their interactions with other 
hadrons are the same at leading order (LO). Heavy quark 
spin symmetry was alread y w idely used in predicting new 
hadronic molecules [7l.ll642~ 



22|. In this Letter, in addition 



to spin symmetry, we will argue that heavy quark flavor 
symmetry is also very useful in the context of hadronic 
molecules, and may be used to verify hadronic molecular 
hypothesis and predict new hadronic molecules. As ex- 
amples, assuming the X(3872) and Zb(10610) to be the 
DD* and BB* molecules, we can predict a series of new 
hadronic molecules, including the Zc(3900) reported very 
recently by the BESIII Collaboration [i^ and later also 
by Belle Q. 

Let us consider the interaction between two heavy 
hadrons forming a bound state. As far as the hadrons 
are not too tightly bound, they will neither overlap nor 
probe the specific details of the interaction binding them 
at short distances. If they contain light quarks, they will 
be able to exchange pions when it is allowed by quantum 
numbers. But the Compton wavelength of the pions is 
far larger than the size of the heavy hadrons or the typ- 
ical length scale associated with the exchange of other 
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mesons, such as the p and uj. That is, there is a scale 
separation between short and long range physics that we 
can exploit to formulate an effective field theory (EFT) 
description of hadronic molecules. 

The situation is analogous to that of the EFT formula- 
tion of the nucleon-nucleon interaction pif , which we use 
as a template for the EFT of heavy hadronic molecules. 
Yet heavy hadrons entail interesting simplifications over 
nucleons. On the one hand heavy quark symmetry heav- 
ily constrains the low-energy interactions among heavy 
hadrons 18 , 2^ . On the other hand pions exchanges are 
20j , in contrast to nuclear physics 
As a consequence, at LO the 



in general perturbative 
where they are not [27 1. 
EFT description only involves energy-independent con- 
tact range interactions 2^, 21 1. 



We are mainly interested in two manifestations of 
heavy quark symmetry: heavy quark spin symmetry 
(HQSS) and heavy flavour symmetry (HFS). Their role 
can be easily illustrated in the heavy meson-antimeson 
system with a series of examples. We begin with HQSS 
as applied to the and Z^, where we assume that they 
are l+~ BB* and B* B* isovector molecules. HQSS im- 
plies that the LO non-relativistic isovector heavy meson- 
antimeson potential is identical in both cases 
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(1) 



where we have indicated the particle channel in the sub- 
script. This explains why the energy shift of the Zi,/Z[^ 
states relative to the BB* / B* B* thresholds is almost the 
same. For a further example we can consider the X(3872) 
- the Xc from now on - as a 1++ DD* molecule. HQSS 
then predicts that the potential in the Xc channel is the 
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identical to that of the 2++ D*D* channel 



) = y^y.42++) 



(2) 



meaning that we can expect the existence of a 2++ HQSS 
partner of the X(3872). Explicit calculations indicate 
that its mass should be in the vicinity of 4012 Me V [2l| . 
Following the previous naming pattern, we will call this 
state the X'^. 

As can be appreciated the exciting feature about heavy 
meson molecules is their high degree of symmetry. This 
is even more evident when we consider HFS. According 
to HFS the interactions involving heavy mesons do not 
depend on the heavy quark flavour. This means that 
the heavy meson-antimeson potential is not able to dis- 
tinguish the D/D* mesons from the B/B* ones. If we 
apply this idea to the Xc, we find 



y^Lg.(i++) = ^B"^-(i++) 



(3) 



and the same is true for the potentials in the X'c Zb 
and Zl channels. The consequence of HFS is that heavy 
meson molecules can appear in flavour multiplets. A res- 
onance in the charm sector might have a counterpart in 
the bottom sector and viceversa. However, there is a 



TABLE I. Various combinations having the same contact term 
as the X(3872) (left) and Zb(10610) (right). Here P and P* 
represent D,B and D*,B*, respectively. 





C'ox 




Ciz 


0(1++) 
0(2++) 
0(2+) 


^.{PP'-P-P) 
P'P* 
D*B* 


1(1+-) 
1(1+-) 
1(1+) 


^,{PP*+P'P) 
P'P* 
D'B' 



catch. The formation of bound states does not only de- 
pend on the strength of the potential, but also on the 
reduced mass of the two-body system. A higher reduced 
mass translates into a stronger binding. If the Xc binds, 
it is more than likely that the Xb - the bottom counter- 
part of the Xc - binds too. Searching for such a state 
may even be regarded as a test of the hadronic molecular 
hypothesis of the X(3872). On the contrary, the shal- 
low nature of the Zb and indicates that their charm 
counterparts are probably unbound. Yet the Zc and Zc 
might survive as virtual states or resonances. As we will 
see, this is indeed the case. 

Now we compute the the expected location of the 
HQSS and HFS partners of the X^ Zb and Z^. For 
that, we notice that at LO the EFT potential is simply 
a contact-range interaction of the type 



(Pl^rip") 

{p\vk°\p') 



Cox 
Ciz , 



(4) 
(5) 



where the subscripts indicate the isospin and whether we 
are considering an X- or Z-like channel (see Table |T]). For 
finding bound state solutions we iterate this potential in 
the Lippmann-Schwinger equation (LSE), where the de- 
tails can be consulted in Ref. [2l[. At this point we find 
it worth commenting that the contact-range potential is 
singular and requires a regularization and renormaliza- 
tion procedure. We employ a standard gaussian regu- 
lator with a cut-off A = 0.5 — 1 GeV, where the cut-off 
window is such that A is bigger than the wave number of 
the states, but at most as large as the breakdown scale 
of the theory 21 1 . 

For the numerical calculations, we work in the isospin 
symmetric limit and use the averaged masses of the 
heavy mesons, which are Mjj = 1867.24 MeV, Md* = 
2008.63 MeV, Mb = 5279.34 MeV and Mb- = 
5325.1 MeV. The value of Cqx is determined from repro- 
ducing the central value of the Particle Data Group aver- 
aged mass of the ^^(3872), 3871. 68±G.17 GeV [28[. The 
resulting value is Cqx = —1.94 fm^ for A = 0.5 GeV and 
—0.79 fm^ for A = 1 GeV [2l[, where the uncertainties 
coming from the error in the mass of the Xc are negligi- 
ble. At this point one may argue that isospin breaking is 
important for the Xc, owing to its closeness to the DqDq 
threshold, but concrete calculations indicate that the ef- 
fect is tiny for spectroscopy [13]. In turn, the value of Ciz 
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TABLE II. Heavy meson-heavy meson combinations having the same contact term as the X(3872) and ^6(10610), and the 
predictions of the pole positions, which are understood to correspond to bound states except if we write "V" in parenthesis for 
denoting a virtual state. The masses are given in units of MeV. 



Vc 




States 


Thresholds 


Masses (A = 0.5 GeV) 


Masses (A = 1 GeV) 


Measurements 


Cox 


0(1++) 




3875.87 


3871.68 (input) 


3871.68 (input) 


3871.68 ±0.17 [28] 






D* D* 


4017.3 


4012+1 


4012+7 


? 




0(1++) 


^,{BB'-B'B) 


10604.4 


105801^ 


10539120 


? 




U(Z ) 


R* R* 


iUDOU.Z 


iUDZD_5 


iU0»4_20 
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U(Z ) 


D* B* 


t OOO. ( 


yooo+4 




? 


Ciz 


1(1+-) 


^(BB'+B'B) 


10604.4 


10602.4 ± 2.0 (input) 


10602.4 ± 2.0 (input) 


10607.2 ±2.0 [5| 
10597 ± 9 [29] 




1(1+-) 


B*B* 


10650.2 


10648.2 ± 2.1 


10648.1 ±2.2 


10652.2 ± 1.5 [5] 
10649 ± 12 [29] 




1(1+-) 


■j^iDD* +D*D) 


3875.87 


38701^0 (V) 


383611? (V) 


3899.0 ±3.6 ±4.9 [23] 
3894.5 ±6.6 ±4.5 [24] 




1(1+-) 


D'D* 


4017.3 


40131;^ (V) 


39831^^ (V) 


? 




1(1 + ) 


D*B* 


7333.7 


7333.51?;^ (V) 


73281^3 (V) 


? 



may be fixed using the Zb(10610) mass. The mass of the 
Zb(10610) measured in the T(nS)'K, hb{nP)'n distribution 
10607.2 ± 2.0 MeV [| is 1.3 a above the B^B* threshold, 
while the value measured in the T(55) — > BB*tt decay 
10597 ± 9 MeV ^ overlaps with the BB* threshold. 
However, these estimations are based on parametrizing 
the Zb and poles as Breit-Wigner. The analysis of 
Ref. [1^ , which overcomes this limitation, suggests that 
the Zb and Zj^ are slightly below threshold and have a 
binding energy of ~ —4.7 MeV and ~ —0.1 MeV respec- 
tively. In line with the estimates of Ref. [s^l , we assume 
the Zb binding energy to be —2.0 ± 2.0 MeV, yielding 



- + 15 
-28 

for A 1 GeV. 



Ciz = -0.75liS fm^ for A = 0.5 GeV and -0.30l? fm^ 



With these values we can make predictions by solving 
the LSE, as previously commented. We summarize our 
results in Table |lTl The uncertainties that are listed cor- 
respond to taking into account that HQSS and HFS are 
not exact, but approximate. We expect a Aqcd/™q de- 
viation of the Cqx and Ciz value from the heavy quark 
limit, which translates into a relative 15% error in the 
charm sector and 5% in the bottom one. Actually, the 
errors are dominated by the uncertainty in the charm sec- 
tor. When we compute the Xb and X^, the relative error 
of Cqx is rather 15% than 5% as its value has been deter- 
mined from the Xc(3872). We remind that the uncertain- 
ties coming from the errors in the mass of the Xc(3872) 
are negligible in comparison. For the states derived from 
the ZbS we sum the Aqcd/ttiq and the binding energy 
errors in quadrature, where the binding error dominates. 
Some of the states - the partners of the Zb/ Z[^ - are not 
bound, but virtual. We indicate this with a "V" . 

Among the predicted states, the 2++ ones can decay 
into two heavy pseudoscalar mesons in a _D wave, which 



would introduce a width of order 0(10 MeV). We refer 
the refined results taking into account the coupled chan- 
nels to a forthcoming work. The predicted mass of the 
D*D* bound state i s hig her than the Xc2(2P) with a mass 
of 3927.2 ±2.6 MeV [28[, and might be searched for in the 
same process as the Xc2(2P), i.e. 77 — > DD. The data 



collected at both Belle and BaBar [3l|, |32[ in that range 
do not have enough statistics for concluding the existence 
of such a state. The most robust prediction would be the 
BB* bound state with I{J^^) = 0(1++), to be called 
Xf,(10580), the analogue of the Xc(3872) in the bottom 
sector. This state should be narrow since the decay into 
the BB is forbidden. It would decay dominantly into a 
bottomium and light mesons. Moreover, the difference 
between the charged and neutral BB* threshold is tiny, 
and completely negligible when compared with the bind- 
ing energy. Therefore, unlike the Xc(3872), whose decays 
exhibit a large isospin breaking, the Xb(10580) would de- 
cay into T{nS)TnnT (n = 1,2) rather than T(n5)7r7r. It 
can also decay into Xb,j{nP) and pious. 

The Zc and Z'^ appear as virtual states, not very far 
away from their respective thresholds. However, the un- 
certainties of the LO calculation are large. In principle 
we are tempted to identify the Zc(3870) state reported in 
Table E] with the Zc(3900) found by the BESIII collab- 
oration. This idea is strengthened when we consider for 
example the possible corrections to the results of Table ITTl 
from coupled channel dynamics. HQSS predicts that the 
Zc and Z'c channels couple, meaning that the location of 
the states will change slightly: the Zc will get a bit of 
extra attraction and the Z'^ repulsion. As a consequence, 
the Zc will increase its mass and eventually bind, while 
the Z'^ Will become deeper and hence unobservable (be- 
ing a virtual state). But without knowing the coupling 



4 



between the and Z^, the estimation of the mass shift 
remains elusive. By assuming the Z^-Z'^ potential to be 
±0.5 C\z we obtain a shift in the position of the states of 
the order of 5 McV for a 1 GeV cut-off. Yet there may be 
other coupled channel corrections available, for instance 
hc{2P)TT or iP{2S)tt, that will require consideration for 
an improved estimation of the mass of the Z^- All this 
indicates that the Zc looks a promising candidate to ex- 
plain the recently observed Zc(3900) resonance, though 
further theoretical effort is still required. 

To summarize, in this Letter we have argued that 
in addition to HQSS, HFS can be used to predict new 
heavy meson molecules. We have also considered the 
uncertainties due to the finite mass of the heavy quarks. 
The predictions might be important in understanding the 
newly observed hadrons in the heavy quarkonium mass 
region in the sense that, if the XYZ states arc hadronic 
molecules, they will probably have heavy flavour partners 
that should be searched for. Note that HFS is a symme- 
try among the coefhcients in the interaction of the La- 
grangians (or, equivalently, the heavy meson potentials), 
and not a symmetry in the binding energies: what mat- 
ters for binding is the potential times the reduced mass 
of the heavy hadrons. Particularly, we studied in detail 
the new states that can be derived from the hypothe- 
sis that the X{3872) and ^^(lOeiO) arc DD* and BB* 
hadronic molecules, respectively 0. Searching for the 
isoscalar 1++ BB* bound state in the T(1S', 2S')7r+7r-7r'' 
channel at hadron colliders or photon-photon collisions 
would provide valuable information on the structure of 
the X{3872). In addition, we find a promising isovector 
1^ DD* virtual state near threshold that might very 
well be identified with the newly discovered Zc(3900) {231 . 
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^ We notice that approaches involving phenomenological (i.e. 
model-dependent) ingredients can lead to other conclusions: 
while the Xi, is usually predicted, the Zc is not [Tsl Issl . |3^ . 
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